The energy spectrum of cosmic rays above 10 15 eV as derived 
from air Cherenkov light measurements in Yakutsk 

A. A. IvanovO S.P. Knurenko, and I.Ye. Sleptsov 
Institute of Cosmophysical Research and Aeronomy, Yakutsk, 677891 Russia 

Abstract 

The Yakutsk array observes the Cherenkov light emitted by UHECR in atmosphere. Recently, 
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an autonomous subarray is added consisting of photomultipliers to measure the showers in the knee 



region. Our aim is to analyze the combined data set in order to derive the cosmic ray spectrum in 
the energy range as wide as possible using the same technique. 

The advantage of the air Cherenkov light measurement is the model independent estimation of 
the EAS primary energy using the total light flux emitted in the atmosphere. A set of the light 
' lateral distributions observed in the extended energy range is presented also. 
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I. INTRODUCTION 



The Cherenkov light emitted in atmosphere by ultrarelativistic electrons of extensive 
air showers (EAS) of cosmic rays (CRs) carries the important information on the primary 
energy and shower development in atmosphere. At the Yakutsk array air Cherenkov light 
measurements are used to estimate the primary energy in the model independent manner, 
and the lateral distribution form of the light on the ground is used to infer the position 
of shower maximum in atmosphere, X max Q]. In this paper we have derived the energy 
spectrum of cosmic rays in the range ~ 1 to 6 x 10 4 PeV using the independent measurement 
technique via air Cherenkov light and compared it with our previous spectrum obtained with 
EAS charged particle measurements as well as the results of other giant arrays. 

II. AIR CHERENKOV LIGHT MEASUREMENTS WITH THE YAKUTSK AR- 
RAY 

During more than 30 years of lifetime the Yakutsk array is detecting muonic, electronic 
and Cherenkov light components of EAS. The actual detector arrangement of the array is 
shown (Fig. Open circles are charged particle detectors given as the background here. 
Filled circles show the Cherenkov light detectors - open photomultiplier tubes (PMTs) of 
176 cm 2 and 530 cm 2 acceptance area, forming a medium subarray. Filled squares indicate 
PMTs of autonomous subarray with independent trigger. This subarray was added in 1995 
with the aim to study air showers in the energy range 10 15 — 10 17 eV via the Cherenkov light 
measurements. 

During ~ 15000 hours of observation in moonless nights there were detected ~ 60000 
showers of energy above 6 x 10 16 eV by the medium array. They were selected using the 
double trigger condition: i) three or more scintillator stations have detected the charged 
particle density above the threshold 1 particle/m 2 ; ii) The Cherenkov light intensities in 
three or more PMTs are greater than 2.4 x 10 5 /4.5 x 10 5 photons/m 2 depending on the 
acceptance area of PMT. 

The autonomous array data consist of ~ 200000 showers with E > 1.2 x 10 15 eV detected 
during ~ 3200 hours of observation. The Cherenkov trigger condition only is used to select 
these showers. 
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FIG. 1: The detector arrangement of the Yakutsk array. A description of the points is given in 
the text. 

For absolute calibration of PMTs the Cherenkov radiation of relativistic particles in 
distilled water was used [lj. The variations in cosmic ray intensity measured with Cherenkov 
detectors itself are used to monitor the atmospheric transparency Q]. 



III. LATERAL DISTRIBUTION OF THE CHERENKOV LIGHT 

We have summarized our lateral distribution function (LDF) measurements of the 
Cherenkov light intensities in Fig. El The data of both subarrays (zenith angle 9 < 30°) 
are parameterized by the light intensity at 150 m from the shower core, Q150, the only core 
distance really present in the shifting range of measurements when the primary energy is 
rising from ~ 10 15 to 10 19 eV. The data are consistent with the previous results of the 
Yakutsk array concerning the Cherenkov light and can be described by the suitable EAS 
model simulation There's seen no abrupt change of LDF parameters, so we choose 
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TABLE I: The fraction of Qtot actually measured (A). 
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rather smooth approximation curve to fit the experimental data in the whole energy range: 

(R 1 + 150)(R 2 + Ry- b 
Q[R) ~ Ql50 (R 1 + R)(R 2 + 150y-»> (1) 

where R l = 60 m; R 2 = 200 m; b = (1.14 ± 0.06) + (0.30 ± 0.02) x lg Q 150 



IV. ESTIMATION OF THE PRIMARY ENERGY OF AIR SHOWERS 

The total flux, Qtot, of the Cherenkov light emitted by the shower in atmosphere is used 
to estimate the primary energy. The detector design of the Yakutsk array is appropriate to 
measure Q to t in the range above about Q150 = 10 7 m~ 2 (Table HJ). 

In each Q150 bin the LDF extrapolation formula Q is used to calculate the total flux. In 
order to connect Q to t with the ionization loss, Ej, of the shower in atmosphere one needs to 
know the atmospheric transparency factor for the Cherenkov light and model predictions of 
the shower development. We refer to |l| where this factor was given on the basis of model 
calculations and transparency measurements in Yakutsk: 

= 2.18 x 10 4 Q tot 
1 0.37+1.1 x 10- 3 X, ' ' ' 



max 



This factor is nearly model-independent in a sense that the model dependence is param- 
eterized by X max . 

The fraction of the primary energy deposited in atmosphere can be evaluated using the 
total number of electrons and muons measured with the Yakutsk array, the attenuation 
length of the shower and the mean energy of muons. The model calculation results are 
used to estimate the energy of EAS nuclear active component and neutrinos. The resulting 
apportioning of the primary energy 10 18 eV is given in Table ITT1 

We have used these values to derive the primary energy, E, of showers with measured 
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FIG. 2: Air Cherenkov light radial distributions. Our experimental data are given by the points. 
The curves are approximations along Equation ^ for given Qiso- 
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TABLE II: Estimation of the primary energy portions gone with the shower components. 



Energy deposit channel 


The portion 




of energy, % 


Ionization loss above sea level 


78 


The total electromagnetic component energy 


87 


Energy transferred to muons and neutrinos 


9 


Energy carried by the nuclear active component 


4 



Ql50 : 

E = (1.07 ± 0.27) x 10- 5 Q[° 50 99±0 - 02) , PeV (3) 

V. THE ENERGY SPECTRUM OF COSMIC RAYS DERIVED FROM AIR 
CHERENKOV LIGHT MEASUREMENTS 

In order to evaluate the intensity of the primary flux we have collected data during the 
observation periods with atmospheric transparency better than 0.65, shower axes within area 
of corresponding subarray, zenith angle 6 < 30°. Acceptance area of the array, S e ff, depends 
on the primary energy and zenith angle. We have calculated S e ff as a function of Q150 
averaged in zenith angle interval (0°,30°) for each subarray (Fig. EJ). Lateral distribution 
fit ^ is used along with instrumental and statistical errors to model the trigger of each 
subarray with 100000 fake showers. In the case of autonomous array the Cherenkov trigger 
is simulated while in the medium case the double trigger for scintillator and PMT signals is 
modeled. The shower data gathered after the array was re-configured were used to work out 
the spectrum. Namely, 1993-2001 for the medium, and 1995-2001 for autonomous subarray. 

The resulting differential all-particle spectrum of cosmic rays is given (Fig. Ta- 
bles 11111 11VJ) together with the Yakutsk array spectrum carried out by the charged particle 
measurements with the trigger-500 m |5]. Spectra are compatible above 2000 PeV within 
errors but are contradictory at lower energies. We have to re-analyze our data near the 
threshold energy of the scintillator trigger-500. Both methods confirm the spectrum irreg- 
ularity near E ~ 10 4 PeV (an 'ankle'), and the autonomous subarray data is compatible 



6 



7 r 

6.5 
6 




logio(Qiso/1 photon m" 2 ) 



FIG. 3: Acceptance area of the Cherenkov subarrays. The solid curve is for autonomous, and the 
dotted one represents the medium subarray. 

with a 'knee' at E ~ 3 PeV. Our next spectrum obtained with trigger-1000 measurements 
of charged particles p] is in the energy range beyond both shown spectra, out of reach of 
the Cherenkov subarrays, therefore is not given here. 

The autonomous subarray spectrum in the range 1 to 100 PeV was given in comparison 



with other relevant results earlier 



Comparison with results of othergiant arrays (Fig. |3J) 
shows a satisfactory agreement of the flux with Akeno Q and AGASA (4J data in the energy 
range covered by the Cherenkov detectors of the Yakutsk array. The data of HiRes-I and 

Q 

HiRes-II from the preprint 7} have the divergence with our spectrum which can be explained, 
presumably, by the difference in the primary energy estimation ~ 30%. Re-evaluated energy 
spectrum of the Haverah Park array 8] in the range 3 x 10 2 PeV to 4 x 10 3 PeV coincides 
with the HiRes-II spectrum and isn't given here. 
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FIG. 4: The cosmic ray flux measured by the air Cherenkov light (squares for autonomous and 
circles for medium subarrays) and charged particle detectors (trigger-500, triangles) of the Yakutsk 
array. Each data point represents the differential flux scaled up by a factor of E 3 . Error bars 
represent statistical errors only. 

VI. CONCLUSIONS 

The independent measurement technique using air Cherenkov light detectors of the 
Yakutsk array enabled us to obtain the cosmic ray energy spectrum in the range ~ 1 to 
6 x 10 4 PeV. Although the PMTs arrangement of the array isn't convenient to study air 
showers above 10 4 PeV, there is an objective field of activity below that border, wide enough. 
Our preliminary 'Cherenkov' spectrum contradicts previous 'scintillator' spectrum around 
the threshold energy of the trigger-500. At the higher energies two spectra agree within 
errors. Both spectra exhibit an 'ankle' feature around 10 4 PeV. The autonomous subarray 
data confirm the 'knee' at E ~ 3 PeV. A comparison of our results with the data of other 
giant EAS arrays shows the similarity of spectra within ~ 30% uncertainty in the primary 
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FIG. 5: Our energy spectrum in comparison with other experiments. The data are from HiRes-I 
(open triangles), HiRes-II (filled triangles), Akeno (filled squares), AGASA (open squares). Our 
results are from the autonomous subarray (filled circles) and medium subarray (open circles). 

energy estimation. 
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TABLE III: CR flux measured with the autonomous subarray. 


Errors represent the statistical 


uncertainties. 




Primary energy bin 


Differential flux J(E) 


\og 10 (E/leV) 


(m _2 sr _1 s _1 eV _1 ) 



15.15 - 


15.25 


8.54 x 10~ 22 


15.25 - 


15.35 


4.63 x 10~ 22 


15.35 - 


15.45 


2.56 x 10~ 22 


15.4 


- 15.5 


1.80 x 10~ 22 


15.5 


- 15.6 


1.05 x 10~ 22 


15.6 


- 15.7 


4.83 x 10~ 23 


15.7 


- 15.8 


2.32 x 10~ 23 


15.8 


- 15.9 


1.68 x 10~ 23 


15.9 


- 16.0 


9.06 x 10~ 24 


16.0 


- 16.1 


2.83 x 10~ 24 


16.1 


- 16.2 


1.60 x 10" 24 


16.2 


- 16.3 


7.27 x 10" 25 


16.3 


- 16.4 


3.33 x 10~ 25 


16.4 


- 16.5 


1.73 x 10~ 25 


16.5 


- 16.6 


(8.77 ±0.42) x 10~ 26 


16.6 


- 16.7 


(4.95 ± 0.23) x 10~ 26 


16.6 


- 16.7 


(4.95 ± 0.23) x 10~ 26 


16.7 


- 16.8 


(2.71 ±0.14) x 10~ 26 


16.8 


- 16.9 


(1.20 ± 0.07) x 10~ 26 


16.9 


- 17.0 


(7.10 ±0.43) x 10- 27 


17.0 


- 17.1 


(3.80 ± 0.25) x 10~ 27 


17.1 


- 17.2 


(8.11 ±0.66) x 10~ 28 


17.2 


- 17.3 


(7.43 ± 0.66) x 10~ 28 


17.3 


- 17.4 


(4.25 ± 0.63) x 10~ 28 


17.4 


- 17.6 


(1.10 ±0.23) x 10~ 28 
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TABLE IV: Cosmic ray flux measured with the medium subarray detectors. 



Primary energy bin Differential flux J(E) 



log w (E/leV) (m^sr^s^eV- 1 ) 



17.1 - 17.3 


9.24 x 10 


17.3 - 17.5 


2.16 x 10 zo 


17.5 - 17.7 


5.91 x 10 


17.7 - 17.9 


1.60 x 10 


17.9 - 18.1 


4.20 x 10 


1 O 1 1 O Q 

18.1 - 18.3 


1 11 , , 10 — .30 

1.11 x 10 


18.3 - 18.5 


(2.62 ± 0.13) x 10 d 


18.5 - 18.7 


(6.38 ± 0.58) x 10~ 32 


18.7 - 18.9 


(9.47 ±2.55) x 10~ 33 


18.9 - 19.1 


(2.93 ± 1.10) x 10" 33 


19.1 - 19.3 


(3.00 ± 1.62) x 10" 33 


19.3 - 19.5 


(6.26 ±4.55) x 10" 34 


19.5 - 19.7 


(3.48 ± 2.90) x 10~ 34 


19.7 - 19.9 


(1.71 ± 1.57) x 10~ 34 
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